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Abstract: Many controlled realizations of chaotic inflation employ pseudo-scalar axions. 
Pseudo-scalars <f> are naturally coupled to gauge fields through cipFF. In the presence 
of this coupling, gauge field quanta are copiously produced by the rolling inflaton. The 
produced gauge quanta, in turn, source inflaton fluctuations via inverse decay. These new 
cosmological perturbations add incoherently with the "vacuum" perturbations, and are 
highly nongaussian. This provides a natural mechanism to generate large nongaussianity 
in single or multi field slow-roll inflation. The resulting phenomenological signatures are 
highly distinctive: large nongaussianity of (nearly) equilateral shape, in addition to de- 
tectably large values of both the scalar spectral tilt and tensor-to-scalar ratio (both being 
typical of large field inflation). The WMAP bound on nongaussianity implies that the 
coupling c of the pseudo-scalar inflaton to any gauge field must be smaller than about 
10 2 M p l . 
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1. Introduction 



Primordial inflation provides a simple mechanism to resolve the conceptual difficulties of the 
standard Big Bang cosmology and has enjoyed great phenomenological success in account- 
ing for the properties of the observed Cosmic Microwave Background (CMB) anisotropics. 
As such, inflation has become the dominant paradigm for the early universe. In spite of 
this success, however, a compelling particle physics realization is still lacking. The key 
obstruction is the requirement of a suitably flat scalar potential, V(<p). Successful inflation 
requires e, \rj\ <C 1 where the slow roll parameters are defined as 

Ml (V'\ 2 9 V" , x 

Here prime denotes derivative with respect to ip, and M p ~ 2.4 • 10 18 GeV is the reduced 
Plank mass. These parameters are notoriously sensitive to Ultra- Violet (UV) physics. For 
example, even generic Planck-suppressed corrections to V((p) may contribute A 77 = 0(1), 
thus spoiling inflation. This UV sensitivity represents a technical fine tuning problem which 
must be addressed in any particle physics model of inflation. 

It is conceivable that dangerous corrections to e, 77 may be absent as a result of fine- 
tuning [|J, or that the requirement of a flat potential can be evaded by invoking somewhat 
exotic effects such as dissipation |2|, |J| , small sound speed Q or higher derivative corrections 
||. However, perhaps the simplest and most cogent way to realize e, |r/| <C 1 in a natural 
way is by assuming that the inflaton <p is a Pseudo-Nambu-Goldstone-Boson (PNGB) 
H |7|, |8|, |9|, [l(], [ll], 12, 13, 14, 15]. In this case the inflaton enjoys a shift symmetry 



(p — > ip + const, which is broken either explicitly or by quantum effects. In the limit 
of exact symmetry we must have e = 77 = 0, thus dangerous contributions to the slow roll 
parameters are controlled by the smallness of symmetry breaking. Moreover, PNGBs like 
the axion are ubiquitous in particle physics: they arise whenever an approximate global 
symmetry is spontaneously broken and are plentiful in string theory compactifications. 

The idea of invoking a PNGB to obtain a natural realization of inflation is more than 
20 years old. The first model, natural inflation ||, exploited the periodic potential 



V(p) = A 4 



1 / v 

1 — cos — 



(1.2) 



which arises from nonperturbative effects and breaks the continuous shift symmetry down 
to a discrete subgroup p — > <p + (27r)/. Unfortunately, this model is compatible with obser- 
vation only when the axion decay constant is / > M p |l6| , a regime that may be impossible 
to realize in a controlled effective field theory because it suggests a global symmetry broken 
above the quantum gravity scale ||. It has also been conjectured that / > M p cannot be 
realized in a controlled limit of string theory 17]. Fortunately, this difficulty can be easily 



evaded, for example by considering two 3 or more [10, [Tl]] axion fields, extra dimensions (8|, 
or by exploiting the non-periodic contributions to V(<p) that arise from wrapping branes on 
suitable cycles [l3|]. Currently, there exist a number of natural, controlled realizations 
of axion inflation with an axion scale / a few orders of magnitude smaller than M p but 
that nevertheless behave effectively as large field inflation models 1C, 11, |l^, 14, 15]. 
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In any axion inflation model, the inflaton is expected to couple to some gauge field 
via interactions of the type 

C ist = ~ V >F^F lu/ (1.3) 

where F^ = d^A u - d u A^ is the field strength and F^ v = ± e ^ Q/3 ' F aj3 is its dual. The 
strength of the interaction is controlled by the decay constant, /, and by the dimensionless 
parameter a. While a is in principle a model dependent quantity, from the perspective of 



effective field theory we generally expect it to be order unity. On the other hand, Ref. [15 



has provided concrete examples (in multi-field or extra-dimensional models) which can 
result in a greater or much greater than one. Thus, for controlled effective field theory 
realizations of axion inflation it is very natural to have a/ f 3> M" 1 , in which case the 
interaction (|l.3|) is much stronger than gravitationally suppressed. In this work, we consider 



in detail the cosmological implications of the interaction (1.3), which has been neglected 
in nearly all previous studies. Our analysis is quite general: the general logic of effective 
field theory requires the inclusion of an interaction (|1.3|) whenever <f> is pseudo-scalar. 



In [15] it was shown that energy dissipation into gauge fields, via the interaction (|1.3j), 
can slow the motion of the inflaton on a very steep potential. A more conservative approach 



was adopted in our recent work [18], where we note that the interaction (|1.3j ) can have a 
profound impact on the phenomenology of the model, even in the conventional slow roll 
regime. The underlying physics is as follows. The motion of the inflaton amplifies the 
fluctuations of the gauge field, SA, which in turn produce inflaton fluctuations via inverse 
decay: 5 A + 5 A — > 5ip. When f /a ^ 10~ 2 M p , which is very natural in models that admit 
a UV completion, this new source of perturbations actually dominates over the usual 
fluctuations from the vacuum. In this regime, all previous studies of axion inflation are 
invalid. Our analysis is phenomenological: we use CMB data to place observational limits 
on the coupling a/f without making any specific assumptions about the microphysical 
origins of the model. We believe that this study should serve as motivation for a case-by- 
case analysis of the allowed values of a in various explicit realizations of axion inflation. 

In this paper, we reconsider the cosmological fluctuations in axion inflation, extend- 
ing significantly our previous work [18|. We are motivated, in part, by the recent surge 
of interest in computing and measuring nongaussian effects in the CMB (see jl{| for a 
recent review). Nongaussian statistics, such as the bispectrum, provide a powerful tool to 
discriminate between the plethora of inflationary models in the literature and may pro- 
vide a valuable window into the detailed physics of the very early universe. However, a 
single decoupled scalar field in slow roll is well known to produce an undetectably small 
signal p0| , |2l| , p^ |. The reason for this is intuitively easy to understand: nongaussianity 
is a measure of the strength of interactions, while in the vanilla scenario the requirement 
of a flat potential typically also constrains interactions to be weak. To evade this no-go 
result and obtain an observably large signal, previous studies have invoked non-standard 
field theories (with small sound speed [23| or higher derivatives ||) or initial conditions 
|23| , ^||, |2f| p6| , potentials with sharp features [27, ^q] , dissipative effects ||, fine-tuned 
inflationary trajectories [29| or post-inflationary effects (such as preheating [30, [H]]). In 
the work fl~8||, on the other hand, the no-go results of [p0|, |2l|, [22| are circumvented in a 
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very novel way: the interaction that gives rise to large nongaussianity, eqn. (1.3), does not 
play any role in the background dynamics and is thus unconstrained by the requirement of 
slow roll. 

In terms of the broader picture we believe that this work - along with flq] and the 



previous studies p2| , 33, 34] - suggests that the conventional lore concerning the difficulty 
of obtaining large nongaussianity may have been excessively conservative. One generically 
expects that inflaton to couple to some fields which do not play any role in driving inflation. 
Such interactions are unavoidable from an effective field theory perspective and are (at least 
to some extent) necessary in order to successfully reheat the universe after inflation. Refs. 



1 18 , |32], |34[ provide explicit examples demonstrating how the consistent inclusion of such 
interactions can radically modify the phenomenology of inflation, via particle production 
effects. This work represents a challenge to the conventional lore that nongaussianity is a 
"smoking gun" signature of non-standard inflationary dynamics by illustrating explicitly 
that perhaps the simplest and best-motivated particle physics models of inflation are already 
constrained by existing observational limits on nongaussianity. 

This paper is organized as follows. In section ^ we provide an overview of the mecha- 
nism and briefly summarize our previous work Jl8| . In section || we compute the correlation 
functions of the curvature fluctuations and tensor (gravity wave, GW) perturbations. Sec- 
tion |] studies the resulting phenomenology. We have written this section in self-contained 
way, so that a reader who is interested only in the observational predictions can skip the 
previous more technical section. In section || we perform a complete computation of the 
perturbations, that includes also the metric perturbations. We show that metric pertur- 
bations can be neglected (at leading order), in the computation of the density and GW 
correlators sourced by the gauge field. In section || we review most of the existing models 
of axion inflation, and briefly discuss the implications of our findings for these models. 
Finally, in section [|, we conclude. 

2. Overview of the Mechanism 

In this section we provide a brief overview of the production of gauge field fluctuations 
and the subsequent inverse decay effects in axion inflation. This section is largely review 



of |J|. 



We consider a simple theory of a PNGB inflaton interacting with a U(l) gauge field 1 



via the interaction (1.3). The action is 



(2.1) 



where R is the Ricci scalar, F^ v = d^A u — d u A^ the field strength, and F^ u = 1 F a p 
its dual, with r/ 0123 = 1. We separate the inflaton into a homogeneous (background) part 
plus its fluctuations 

ip = (f)(t) + 6(p(t, x) (2.2) 



1 The generalization to non-Abelian gauge groups is straightforward. 
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We leave the potential V(cp) arbitrary, except to assume that it is sufficiently flat to sup- 
port the required amount of inflation (iV e > 60). We assume a spatially flat Friedmann- 
Robertson- Walker (FRW) space-time with metric 

ds 2 = g^dx^dx" = -dt 2 + a 2 {t) dx ■ dx (2.3) 
= a 2 (r) [-dr 2 + dx-dx] (2.4) 

where on the second line we have introduced conformal time, r, related to cosmic time as 
adr = dt. Derivatives with respect to cosmic time are denoted as dtf = / and with respect 
to conformal time as d T f = f. The Hubble rate H = a/a has conformal time analogue 
H = a! I a. 

We are first interested in the gauge quanta which are produced by the homogeneous 
rolling inflaton <p(r). To this end, we can ignore the inflaton and metric perturbations 
in the equations of motion of the gauge field (see section [5] for the complete treatment). 
Extremizing the action with respect to Aq, and choosing the Coulomb gauge Aq = 0, then 
gives • Aj = 0, from which we set V • A = 0. The equations of motion for A then read 

A" -V 2 A-j(j)'V x A = (2.5) 



As we discuss in subsection 2.1, this equation describes the production of the quanta of 
the gauge fields that results from the motion of the inflaton. 

The produced gauge quanta have two key effects: they backreact on the homogeneous 



background dynamics (see subsection |2.2| ) and also source inflaton perturbations (see sub- 



section 2J3). Both effects are governed by the equation of motion of the inflaton, and the 



00 Einstein equation, which read, respectively 

if" + 2-Hif' -VV + a 2 ^ = a 2 -E • B 

dip f 



H 2 - ' 



3M p 2 



i^ + i(v,) 2 + ^4(m 



(2.6) 



In these equations we have retained the spatial dependence of <p (due to the inflaton 
perturbations), and we have introduced the physical "electric" and "magnetic" fields 2 

B = ^Vxl, E = -\l' (2.7) 
a z a A 

2.1 Production of Gauge Field Fluctuations 

During inflation, the motion of the inflaton leads to an instability for the fluctuations 
of the gauge field. To see this effect, we start from the equation of motion for A^ in 



the background of the homogeneous inflaton 4>(t), eq. ( |2.5D above. We decompose the 
q-number field A(t, x) as 

d?k 



j-^ [6 A (k)a A (k)yl A (r,k) e 4k - x + h.c. 



2i 



2 We do not assume that necessarily corresponds to the Standard Model electro-magnetic gauge 
potential. 
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where "h.c." denotes the Hermitian conjugate of the preceding term, the annihilation/creation 
operators obey 

a A (k),4,(k')] = W (3) (k - k') (2.9) 
Here e\ are circular polarization vectors satisfying k ■ e± ( k ) = 0, k x e± ( k ) = ^fike± Ik), 



(+ ( —kj = e± [kj , and normalized according to e\ ykj ■ e\> [kj = 5\\i . 

Inserting the decomposition (|2.8|) into eq. (|2.5D results in the equation of motion 



+ k 2 ± 



2k£ 



A±(r,k) = 0, e = 



no 



2fH 



(2.10) 



for the c- number mode functions A±. During inflation the parameter £ may be treated as 
constant, as its time variation is subleading in a slow roll expansion. 

From equation ( [2.10 ) we see that one of the polarizations of A\ experiences a tachyonic 
instability for k/(aH) < 2£. Without loss of generality, we assume that cj> > during 
inflation, so that the mode exhibiting the instability is A + . In appendix A we review the 
solutions of ( [2.10D and show that the growth of fluctuations is well described by 



A+(T,k) 



1 



k 



2k \2&H 



(2.11) 



in the interval (8£) _1 <, k/(aH) < 2£ [^g] of phase space that accounts for most of the 
power in the produced gauge fluctuations. The phase space of growing modes is non- 
vanishing for £ j> 0(1), which we assume throughout. Notice the exponential enhancement 
e n ^ in the solution (|2.11 ), which arises due to tachyonic instability, and reflects significant 
nonperturbative gauge particle production in the regime £ ^ 1. On the other hand, the 
production of gauge field fluctuations is uninterestingly small for £ < 1. Note also that the 
other polarization state, A_(t, k), is not produced and can therefore be ignored. 

We have thus seen that the motion of the homogeneous inflaton <j){t) leads to produc- 
tion of gauge field quanta SA^. There are two key physical effects associated with the 
interactions of these produced quanta with the inflaton. The first effect is the backreaction 
of the produced quanta on the homogeneous dynamics of <fi(t), a(t). In the next subsec- 
tion we study the conditions under which backreaction effects are negligible. The second 
key physical effect is the production of inflaton fluctuations via inverse decay; this is the 
subject of subsection 

2.2 Backreaction Effects 

Backreaction effects can be accounted for using the mean of the field equations fl2.6|) : 

(2.12) 



a 



+ 3H<P + V'((/)) = j{E-B) 



3H 2 



1 

Ml 



l -^ + V{ (t) ) + l -{E 2 + B 2 ) 



(2.13) 



where we have switched to physical time. The expectation values appearing in 

(E.12IJ2.13D 

encode the backreaction of the produced gauge quanta on the homogeneous dynamics of 
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(p(t), a(t). From (|2j|) and we have 



(E-B) 

-(E 2 + B 2 ) 

2 \ / 



1 



4vr 2 a 4 
1 



dA; A; 3 — I A 



4?r 2 a 4 



dA; A; 2 [| 



A' | 2 + A: 2 |A. 



(2.14) 



Since we are studying the backreaction of the produced quanta, we should disregard the 
modes that do not experience the growth discussed in the previous subsection, namely 
all modes A_ and the large momentum (k > 2^aH) modes of A+. As can be seen from 
( 2.10 ), modes of A + with k/(aH) S> 2£ remain in their vacuum state and do not experience 
any tachyonic instability. Such modes contribute to the vacuum energy of the U(l) field 
(in eqn. (2.13); an analogous consideration applies to (E ■ B) in eqn. ( [Hp ) that wc 
assume is canceled by a bare vacuum energy, as is customary in QFT. (In short, we have 
nothing to add to the cosmological constant problem.) This prescription provides a UV 
cutoff k/ (aH) < 2£ in the integrals ( p. 14 ); see the paragraph after eq. fl2.10p . From a 
direct inspection of the solutions of ( p. 10 ), we verified that the integrals (2.14) converge 
in the infrared k — > region, and that they receive their support almost entirely from the 
region (8£) -1 ^ k/{aH) ;$ 2£ in which ( |2.11 ) is a very good approximation of the exact 
solution. Therefore, ( p. 14] ) can be evaluated by using the expression ( |2.11| ) for A+, and 
by integrating only over this region of momenta. In fact, we verified that the momentum 
interval can be extended from to oo, since the expression ( 2.11 ) rapidly decreases outside 
the (8£) -1 < k/(aH) < 2£ interval, and the contribution of the "outer" regions to (PH) 
can be neglected. Proceeding in this way allows for an analytic result: 

(E-B)^-2A-W- 4 ^e 2 < , ( — + — ) ~ 1.4 • 1(T 4 ^ e 2 ^ (2.15) 

£ 4 2 £ d 

This procedure follows Ref. fjD|, where these expressions were first derived. We expect 
that any sensible renormalization prescription will yield results in agreement with fl2.15|) . 

From ( |2.12 .2.13) we can distinguish two distinct kinds of backreaction effects. First, 
the gauge field fluctuations are produced at the expense of the kinetic energy of cp(t). 
This contributes a new source of dissipation into the homogeneous Klein-Gordon equation 
( 2.1 2|) . In order to trust the usual slow-roll inflationary solution, we require that |3iJ<^| = 
| - V'{<j))\ > \ f(E-B)\. Using ( ^15|) , this condition reads 

^<13e 3/2 e^ (2.16) 

2tt|0| 

The second kind of backreaction effect arises because the energy density in produced 
gauge field fluctuations contributes to the Friedmann equation ( |2.13| ). To ensure that the 
expansion of the universe is dominated by the potential energy of the inflaton we require 
3M 2 H 2 = V » \{E 2 + B 2 ). Once eq. ( |2.15| ) is taken into account, this condition reads 



H 



< 146^ 3/2 e"^ 



(2.17) 



Taken together, the constraints ( |2.1(: , 2~T7| ) ensure that the unstable growth of gauge field 
fluctuations does not modify the usual homogeneous inflationary dynamics. 
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2.3 Inverse Decay Effects and Inflaton Perturbations 



Even when (2.16,2.17) are satisfied, the coupling ipFF may still have a profound impact 
on the cosmological fluctuations in the model (p.l|). The perturbations of the inflaton are 
described by the equation 

d 2 



g T 2 ^ Q T 



V 2 + a 2 m 2 



> a 



6<p(r,x) =a z j (E-B-(E-B 



(2.18) 



where m 2 = V" . The solution of ( 2.18| ) splits into two parts: the solution of the homoge- 
neous equation and the particular solution which is due to the source. Schematically, we 
have 

<fy>(r,x) = ^ vac (r,x) + 5v9 inv . deca y(r,x) (2.19) 
V v ' > „ ' 

homogeneous particular 

The homogeneous solution corresponds, physically, to the usual vacuum fluctuations from 
inflation. The particular solution, on the other hand, can be interpreted as arising due to 
inverse decay processes 5A + 5A —> 5ip. This new source of inflaton fluctuations contributes 
directly to the observable curvature perturbation on uniform density hyper surf aces, owing 
to the relation ( ~ — ^r5(p. It was shown in |l8| that the inverse decay contribution to 
the cosmological fluctuations may actually dominate over the usual vacuum fluctuations 
in the regime / W~ 2 M p . This radically modifies the phenomenology of axion inflation. 
In particular, the inverse decay contribution to the primordial cosmological fluctuations 
is highly nongaussian; this is evident already from ( p. 18 ) since the particular solution 
^inv.decay is bilinear in the gaussian field 5A^. 



To proceed with the computation, we decompose 



Sf(T, x) 



d 3 k Q k (r) 



ikx 



air 



(2.20) 



(2tt) 3 / 2 

Eq. ( 2.18j ) then results in (we note that the last term of ( |2.18 ) has no effect on the mode 
functions with momentum different from zero) 3 



d 2 + k 2 + a 2 m 2 - — 



Qk(r) = J k (r) 



Ju(r) 



a 3 (r) 



Q: 



d 3 k 



fj (2vr)3/ 2 



e~ ik x E ■ B 



We separate the mode functions of the two terms in eq. ( |2.19| ) as 

QUr) = Qr(r) + Qr ACCay (r) 

and we discuss each contribution separately. 
The homogeneous term is expanded as 

CT = ^W + ft t (-kK(r) 



(2.21) 
(2.22) 

(2.23) 



(2.24) 



3 In Section g we present the complete computation, including also scalar metric perturbations. We show 
explicitly that equation (2.21) still holds, with only the addition of subdominant (Planck-suppressed) terms 
in the source Jk- 
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The inflaton ladder operators obey 



6(k),6 t (k')l = <5 (3) (k-k') 



(2.25) 



and commute with the ladder operators of the gauge field 



[6(k),a A (k')] = 6(k),4(k'^ 



(2.26) 



The properly normalized homogeneous solutions of ( 2.21 ) are given by the well-known 
result _ 

Mr) = i^V^H^i-kr), + 0(e,r,) (2.27) 

where we have chosen the (arbitrary) phase so that l -fk{ T ) is rea l i n the limit —kr — > 0. 
The vacuum modes ( |2.27| ) are employed in the retarded Green function associated with 



G k (r,r') = i@(r - r') [^(r)^(r') - ^(r)^(r')] 



(2.28) 



which obeys 



d 2 + k 2 + a 2 m 2 - ^ 



G k (r,T r ) = 8(T-T'). 



Using the Green function ( 2.28 ) the particular solution of ( |2.21 ) takes the form 



/oinv. decay / 



dT'G k (T,T')MT') 



(2.29) 



where the source term was defined in eq. ( [2.22 ). We note that this particular solution is sta- 
tistically independent of the homogeneous solution (|2.24| ). In fact, the particular solution 
can be expanded in terms of the annihilation/creation operators a\(]s.),a\(k) associated 
with the gauge field, while the homogeneous solution can be expanded in terms of the an- 
nihilation/creation operators 6(k),6^(k) associated with the inflaton vacuum fluctuations. 
As we pointed out, these two sets of operators commute with one another. 

We are now in a position to compute the correlation functions for the perturbations 
5ip. We present this computation in section ||. In section |] we discuss the resulting 
phenomenology. Finally, in section ||| we show that these result are valid also when metric 
perturbations are also consistently taken into account. 



3. Correlation Functions 



In this section we compute the main phenomenological signatures of the model (2.1). Specif- 
ically: in subsection |3.1| we compute the two point correlation function of the density per- 
turbation C = — "T <Vi i n subsection 3.2 we present the corresponding power spectrum; in 



subsection [3^ we compute the three point correlation function of (; in subsection 3.4 we 
present the corresponding bispectrum; in subsection |3.5| we finally summarize the compu- 
tation of the power spectrum of the gravity waves (GW) modes. In this section we neglect 
scalar metric perturbations for simplicity, however, in section |5| we show explicitly that 
their consistent inclusion does not modify our results. 
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3.1 Two-Point Correlation Function 

We start from the relation £ (r, x) = — Sip (r, x) between the inflaton perturbations, and 
the curvature perturbation on uniform density hypersurfaces. 4 We decompose the latter 
as 



C(r,x) 



d 3 k 



Ck(r)e 



kx 



(3.1) 



(2tt) 3 /2 

so that Ck = — -t As we discussed in subsection |2.3j the inflaton perturbations comprise 

<p 1 1 

of two terms, one being the vacuum fluctuations, and one the fluctuations sourced by the 
gauge quanta; since these two terms are statistically independent, we have 

H 2 



(CkCk 



,2^2 



(Qi ac Qi> c ) + (Q 



inv. decay /^inv. decay 



(3.2) 



The contribution from the vacuum modes is standard. Using eqs. ( 2.24| ), ( |2.25| ), and ( 2.27] ), 
one obtains the well-known result 

n s -l -, 

(^ k n- — (— ) F ^ (3) (k + k') 



2^2 \^hj 



fc 3 



aH 



(3.3) 



in the late time / large scales limit, -S;r < 1. In the second line of ( |3.3| ) we have introduced 

n s = l + 3-2v = l + 0(e,r]) (3.4) 



pl/2 



l + 3-2is = l + 0(e,r]) 
H 2 



2vr|0| 



(3.5) 



We now compute the 2-point correlator of the particular solution ( |2.29 ) 

H 2 f . . „G k (r,r') G k , (r, r" 



/Am. decay / \ Ainv. decay / \\ 

(C k ( r ) C k ' ( T )) 



dr dr ' 



(Jk(r') J v {t")) (3.6) 



a(r) a(r) 

(In the following, we temporarily omit the superscript 'inv. decay' for notational conve- 
nience.) In evaluating this expression, we can make some approximations for the mode 
functions tp k (t) and ip k {r') appearing in the Green function (|2.28 ). Since we are in- 
terested in the power spectrum of modes well outside the horizon we can use the small 
argument limit — kr <C 1 for the modes ( 2.27| ) entering in the Green function ( 2,28j ), 

a{r)H 1 



(t) 



-krY 



for — kr <C 1 



(3.7) 



^2 k 3 / 2 

Notice that we disregard the small slow roll corrections to the amplitude, but we retain 
them in the momentum dependence (since this controls the departure of the spectrum from 
scale invariance). This gives 



(Ck(r)Ck' (r)> 



2H 4 (-kr) 



n s — 1 



A: 3 



dr' dr" Im [ip k (r')] Im [ip k (t")] (J k (r) J k (r")) 

(3.8) 



4 See section ^| for more details on this relation. 
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In this relation, we have already used the fact that the correlator is nonvanishing only for 



|k| = |k'|. We stress that, in ( |3.8| ) we do not employ the approximation (3.7) for the modes 
( Pk( T ') which appear under the integral, since r' must be integrated over. 

To evaluate ( |3.8[ ), we require the source correlator (Jk (V) Jk (t")). Explicit evaluation 
gives 



a 2 ^ 3 )(k + k') f d 3 q 



(r')a(r")J (27T) 



1 + 



q| 2 — q • k 
|q| |k- q| . 



n 2 



:A[r', |q|,|q-k|] A* [r" , |q|, |q-k|] (3.9) 



where 



A[r', |q|, |q-k|] =|q|^ + (r',|q-k|) A + (r', |q|) + |q - kL4' + (r', |q|) A + (r', |q-k|) 

(3.10) 

and where the relations 



|e+(q)-e+(k-q)| 

have been used. 

Putting all together, we have 



1 



1 + 



q| 2 — q • k 
|q| |k- q| 



6 A (-k) = el(k) 



(3.11) 



^!(z^,(3 )(k+k 

4/202 k 3 V 



(Ck(r)Ck' (r)) 



d 3 r 



(2^ 



1 + 



|q| 2 — q • k 
|q| |k — q| 



dr' (-/) Im [p fc (r')] A [r' |q|, |q - k|] (3.12 



Following similar arguments to those made in subsection 2^, we can use the approx- 
imation ([A-3j ) for the gauge field mode functions (we discuss this step more in details at 
the end of this subsection). This leads to 



(Ck(r)Ck' (r)> 



a 2 H & ^< (-jfer) 



n 3 -l 



2 8 ^3 p^2 k 3 



8^ (k + k' 



x / d 3 q* 



1 + 



q* • k 



I q* | | k q* I 



1 1/2 



q*-&| 1/2 Iq^'^ + lq*-^ 



1/2 



1 1/2 



xX 2 



(3.13) 



where the integration variable q* = q/|k| is dimensionless, and where 



l[z] 



— kr 



dxx 3/2 Re ( 



(3.14) 



(notice that x = —kr). As we are interested only in super horizon modes, — kr <C 1, we 
can set to zero the lower extreme of integration of I. It is then manifest that we can set 
v = 3/2 in the argument of the Hankel function, since the slow roll corrections appearing 
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there only modify (in a negligible amount) the amplitude of the correlator, but not its scale 
dependence. This leads to 



poo 

X (z) ~ / dx (sinx — xcosx) e~ z ^ 
Jo 

For future convenience, we rewrite the correlator ( 3.13| ) as 



(c inv.dccay (r) ^nv. decay (r)) _ ^ (^^-l^^^) (k + ^ 



2vr_ 2 

A; 3 

where V was defined in ( |3.5| ) and the dimensionless function /2(C) is 

2 



(3.15) 



(3.16) 



xX 2 



1 + 



q* • A; 



la I 1 / 2 | a 
1 4* I 1 4* 



A- 1 



1/2 



|q*l 1/2 + |q* 



k\ 1 ' 2 



2^ v1qT+ vlq* 



(3.17) 



In general, the function /2(C) needs to be evaluated numerically. However, a simplifi- 
cation is achieved when the argument of I is much greater than one. 



l[z] 



00 x 3 



-Zy/X 



3360 



for z 3> 1 



(3.18) 



As \/|q*| + y |q* — k\ > 1, this approximation is certainly appropriate at large £. One is 
then left with a two dimensional integral that numerically evaluates to 

7 5-1 D~ 5 

mo = ^ » e » 1 ( 3 - 19 ) 

The degree of accuracy of this approximation can be seen in Figure |l]. It is also useful to 
have a fit for fi in the range 2 < £ < 3, as this is the most relevant one for phenomenology 
(as we discussed in section |j). The best monomial fit to j% in this range is 



mo 



3 • 10- 



5.1 



2 < i < 3 



(3.20) 



We conclude this subsection with the justification of the use of (2.11) into the integral 
( 3.12| ). We know that the expressions ( 2 . 1 1[) are accurate approximations to the modes 
A+ (r, k) in the range l/8£ < — r k < 2£. As discussed in subsection 2.2, for any given 
mode k, this interval corresponds to the times for which the amplification of the gauge field 
is maximal. The key point is to ensure that there is a common region in the integration 
space of ( |3.12 ) for which both the functions A + that enter in the expression of A can be 
approximated by ( |2,11|) . This requires that 



or, equivalently, 



— <-|q|r'<2f and — < -|q - k| r' < 2f 



—: < |q*| x < 2£ and — < |q* — k\x < 2£ 



(3.21) 
(3.22) 
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Figure 1: Exact evaluation and large argument approximation of the function ji- 



(x = — kr 1 is the integration variable in ( 3.14| ), and q* = q/|k|). For the super-horizon 
modes that are relevant for phenomenology, we know that x extends from to oo. There- 
fore, for any value of q*, there are always values of the rescaled time x for which either 



A+ (t\ \qJ) or A+ W 



q*|j or A + It', |q* — k\) is maximal, and approximated by (|A-3|) . Only when 
q* — k\, both conditions are valid for the same values of x, and the the approxima- 



tions ( 2.1 1| ) can be used in the whole integrand of ( |2.2| ). Therefore, our result is correct 
only if the integrand of (2.2) is strongly peaked at |q*| = 0(1). We have verified with 
direct inspection that this is indeed the case (We have verified this claim also using the 
representation (|A-2|) of the gauge field modes which is valid arbitrarily deep in the IR; see 
Appendix C.) 

There is a clear physical reason why the integrand is strongly peaked in the region of 
phase space where the conditions ( 3.21| ) are satisfied. Notice that, for the values z = (1) 
which are relevant for the present computation, the expression fl3,15| ) has most of its support 
at x = (1). This means that the "imprint" of the fluctuations from the gauge modes 
occurs when the wavelength of the fluctuation is of the order of the horizon scale (x of 
order one), and it is caused by the inverse decay of gauge field modes, whose wavelength 
is also of the order of the horizon scale (|q*| of order one). This is a natural outcome 
of causality/locality. Identical considerations apply in the computation of the three point 
function that we perform in subsection [3.3j. 



3.2 The Primordial Power Spectrum 

The two point correlation function in momentum space is related to the power spectrum 
by the standard expression 



(CkCk')=PcW^ (3) (k + kO 



(3.23) 
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As we have seen in the previous subsection, £k = Ck &c + ^ nv deca y j anc j ^ wo terms 
are uncor related. From ( |3.3D and (|3.16| ) we get the power spectrum at late times 



PAk)=V 



n s -l 



1 + ^/2(0 e 4w « 



(3.24) 



where V 1/2 = 

3.3 Three-Point Correlation Function 



In this subsection we calculate the three point function of £, 



inv. decay 



tt3 p 

(Cki (r) Ck 2 (r) Ck 3 (t)> = --yj / dr 1 dT 2 dT i 



G kl (t, n) Gfc 2 (t, r 2 ) g fc3 (r, r 3 ) 

a(r) a(r) a(r) 

x (J kl (n) «/k 2 (r 2 ) J k3 (r 3 )) (3.25) 



Proceeding as in subsection to, we arrive to 



(Cki Ck 2 Ck; 



a 



; H 9 



S <3, (kl+k2+k3) /_^ 



P4> s kfklkl 

x [e (q) • e (kx - q)] [e (q - k x ) • e*(-q - k 3 )] [?(q + k 3 ) • e (-q)] 
,0 

xITj / dn [kiTi cos (fciTj) - sin {kin)] 
J — 00 

X.4 [n, |q|, |ki - q|] A [r 2 , |k x - q|, |k 3 + q|] A [r 3 , |k 3 + q|, |q|] (3.26) 

where A was defined in ( |3.10| ). In this expression we have used the fact that the mode 
functions ^4+ are real (which is true both in the approximations (|A-2j) and (2.11)), and we 
have disregarded any (mild) scale dependence. 5 

The correlator depends on the size and the shape of the triangle formed by the vectors 
ki. We denote 

|ki| = A; , I lt2 1 = x 2 k , |k 3 |=X3&; (3.27) 
and, for future convenience, we parametrize the correlator as 



(Ck! Ck 2 Ck 3 



CY 



>H 9 e 6 < 5(3) (ki + k 2 + k 3 ) 1 + x\ + 



,1 ■> 



80(2vr) 7 / 2 / 3 3 e 



fc 6 



3 3 
x 2 x 3 



h (£, x 2 , x 3 ) 



10 



(2tt) 5/2 P 3 e 6 ^ 



, 5( 3 ) (kx + k 2 + k 3 ) 1 + x\ + . 



k 6 



h{^x 2 ,x 3 ) (3.28) 



2 ^3 



We proceed from eq. ( |3.26| ) as we did in subsection 3J.. Using the last expression, eq. 
( 3.10 ), and eqs. ( |A-3| ), we arrive to 



h (£; x 2 , x 3 ) = — 



e 



3vr x 2 x 3 [l + x% + 



J d 3 q* |ct* 1 1/2 |*i - q*| 1/2 |q* + X 3 k 3 \ 1/2 



e (q*) • e(ki ~ q*)] [|q*| 1/2 + \k - q*| 1/2 ] X (|q*| 1/2 + \k - q*^ 



3 Specifically, we have set v — 3/2 in the mode functions entering in the Green function (2.28)). 
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e q* - h) • e -q* - x 3 k 



\k - q*| 1/2 + |q* +X 3 k 3 \ 1/2 



1 



'(q* 



2 A/^(l A; i-q*l 1/2 + |q* + ;E 3 hi 1 '' 2 



+ x 3 k 3 ) ■ e(-q* 



|q* + x 3 fc 3 | 1/2 + |q*l 1/2 

2? 



2t/ — (|q*+x 3 & 3 | 1/2 + |qJ 1/2 

2-3 



(3.29) 



To continue, we can set 
k = (1, 0, 0) 

^3^3 = - 2 (l ~ x 2 + ^3> V- (1 - a: 2 + 3^3 ) (1 + x 2 - x 3 ) (1 - x 2 - x 3 ) (1 + x 2 + x 3 ), 

(3.30) 

and we note that, to a generic vector k = k (s'm9 cos</>, sin 9 sin <fi, cos 9) , corresponds the 
polarization operator 

1 



e + (k) = —= (cos # cos < 
\/2 



z sin q>, cos sin + i cos < 



sin^) 



(3.31) 



(it is immediate to verify that this expression satisfies all the properties listed after eq. 

In the reminder of this subsection we evaluate f 3 for the equilateral configuration 
X2 = x 3 = 1. In the large £ limit, we can use the analytic approximation ( |3.18| ) for Z, and 
perform the momentum integral numerically. We obtain 

os i n— 7 

/ 3 (£;i, i) = — - 9 — , £»1 (3.32) 

The degree of accuracy of this approximation can be seen in Figure ^. It is also useful to 
have a fit for f 3 in the range 2 < £ < 3, as this is the most relevant one for phenomenology 
(as we discuss in section 0). The best monomial fit to f 3 in this range is 



7.4 • 1CT 8 
/ 3 (£ 1, 1) = , 2 < £ < 3 



(3.33) 



3.4 The Bispectrum and Nonlinearity Parameter 

A popular parametrization of nongaussianity is the nonlinearity parameter Jnl, introduced 
by assuming that the curvature perturbation may be expanded as 



C(x)=C 9 (x) + -/ NL [Cl (x) - (C 2 (x))] 



(3.34) 



where Cg( x ) is a gaussian random field (see [^] for a careful discussion of sign conventions). 
Both ( and Q g may be decomposed as in ([O]) so that the relation between the q-modes of 
the Fourier decomposition is 

d 3 p 



3 f d*p 

Ck = C 3 ,k + g JNL y ^3/2 C 9 ,k C 9 ,k-p 



(3.35) 
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Figure 2: Exact evaluation and large argument approximation of the function f% (£; 1, 1) 



By definition, the three point correlator of C vanishes. However, due to the quadratic term 



in ( 3.34 ), the three point correlator of C, is nonvanishing, and can be expressed through a 
sum of two point correlators of ( g . One finds 6 



(Ck! Ck 2 Ck 3 



10 



(2k) 



5/2 



/ NL iH*0 2 <5 (3) (ki+k 2 + k 3 ) §^ 



(3.36) 



where the power spectrum was defined in ( |3.23j ). To obtain this expression, recall that the 
ladder operators are normalized according to ( p. 25 ), and that the power spectrum P (k) 
is related to the two point function as in ( |3.23| ). One should also identify the two point 
function of £ with that of Q g (as the difference in subleading in a perturbative expansion), 
and disregard the the mild scale dependence of the power spectrum. 

By comparison to ( pep , one may define an "effective" (momentum dependent) non- 
linearity parameter, even when the intrinsic nongaussianity is not of the local form (3.34). 



For axion inflation, using the parametrization (3.28) of the three point correlator, we can 
write 

/ 3 (C;x 2 , x 3 )P 3 e 6 ^ 



Znl(£;z2,x 3 ) 



p c (ky 



(3.37) 



where we recall that V 1 ! 2 



2-ttU 



3.5 Power spectrum of the Tensor Modes 

The produced gauge quanta also source tensor metric perturbations (gravity waves). The 



total GW power spectrum was first given in [18]. Ref. |56[ then pointed out that the chiral 
nature of the GW produced by the gauge modes can be probed through the resulting 



6 Note that the factors of (27r) n/2 differ from ||]. This stems fr om a different convention for the nor- 
malization of the Fourier transform (3.1). 
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nonvanishing (BE) and (BT) correlators of the CMB, as studied in |37], |38|. For the 
present model, a positive detection of parity violation would only be possible in a cosmic 
variance limited experiment, and for a limited portion of the parameter space [36]. In 
particular, one needs to be in a regime in which the GW production from the gauge modes 
dominates over that from the vacuum. In the minimal version of the model studied here, 
this region is ruled out by the nongaussianity limit that we discuss in the next section [18]. 
Ref. [36] circumvented this problem by considering the presence of many (AT <L 10 3 ) gauge 
fields, or by the use of a curvaton. 

In the short paper [18] we only reported the final result, reserving the present work for 
the presentation of the details of the computation. As in the meantime this computation 
has been presented in details in ref. [36|, we only provide a quick summary here. The 
tensor modes enter in the spatial components of the metric as gij = a 2 (5ij + hij), with 
ha = dihij = 0. From the Einstein equations, one then finds 



1 

2^ 



d 2 T + — d T - dl) h 



Ml 



-EjEi — I >, I > 



(3.38) 



where TT denotes the transverse and traceless projection of the spatial components of the 
energy-momentum tensor of the gauge field. The computation of the GW production is 
performed analogously to that of the density perturbations that we have presented in details 
in the previous subsections. As for density perturbations, the GW modes produced by the 
gauge quanta are uncorrelated with those from the vacuum, and the two contributions 
add up incoherently in the power spectrum. The two GW helicities are obtained from the 



projectors IL^-^/^ (k) = (k) (k). One finds the two power spectra 



P, 



H 2 



L/R 



7r 2 M p 2 



1 + 



2H 



fh,L/R (£) e 



4tt£ 



(3.39) 



where 



fh,L/R 



d 3 q* 
(2^? 



( 1 ± cos e) 2 1 - q* cos 9 ± y/l - 2q* cos 9 + q 2 



16(1 



dxxfx [sin x — x cos x] 



2C 

x 



q* 



2q* cos 9 + q 2 



Vl*-f.l] 



(3.40) 



and where and 9 are, respectively, the magnitude of the (dimensionless) integration 
momentum q*, and the angle between this vector and the momentum k of the mode. 

In eq. ( |3.39j ), tit = —2e. We note that the contribution to the spectrum of the modes 
from the vacuum, and of the modes sourced from the gauge field have the same scale 
dependence. The reason for this is that the scale dependence of the second term originates 
from the homogeneous solutions of ( p.38 ) - which are the vacuum solutions - employed in 
the Green function. For the same reason, also the two terms in the scalar power spectrum 
have identical scale dependence, see subsection 3.1. 

At large £, the integral over x can be performed as in the previous subsections, and 



one finds [361 



fh, 



4.3 • 10" 



f, 



h,R 



9.2 • 10 



-10 



£ » 1 



(3.41) 
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The degree of accuracy of this approximation can be seen in Figure [|. It is also useful to 
have a fit for fh t L i n the range 2 < £ < 3, as this is the most relevant one for phenomenology 
(as we discussed in section [|). The best monomial fit to fh,L in this range is 



2.6 • 10" 



5.7 



, 2<£<3 



(3.42) 
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Figure 3: Exact evaluation and large argument approximation of the function fh.L ((,)■ 

From ( |3.41| ) we see that, for the purpose of computing the tensor-to-scalar ratio, one 
can disregard the right helicity G W modes produced by the gauge fields. One then obtains 
the result 

o 2H 2 fk^ nT 

PgW = Ph,L + Ph,R 



it 2 Ml \k 



(3.43) 



first reported in 



4. Phenomenology of Axion Inflation 

From the observational perspective, the key quantities which characterize any model of 
inflation are the spectrum of scalar and tensor perturbations, and Pgw-, along with the 
bispectrum of scalar perturbations that encodes the leading departures from gaussian 
statistics. The explicit computations of these quantities was presented in Section ||. In this 
Section, we study the resulting observational signatures. 



4.1 COBE Normalization and Spectral Tilt 



In subsections 34 and 3^ we found that two uncorrelated terms contribute to the power 
spectrum in axion inflation. These are the usual fluctuations generated from the vacuum, 
along with the modes produced by the inverse decay of the gauge quanta excited by the 
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motion of the inflaton. Taking into account both contributions results in a power spectrum 
of the form 



V — 

k 



n s -l 



l + Vf 2 (0 e 4 ^ 



where n s 
defined 



1 + 27/ 



P c (k) 

6e is the spectral index, the pivot is ko 
p l/2 H 2 



2vr|0| 



(4.1) 

0.002 Mpc" 1 , and we have 
(4.2) 



It is worth noting that both terms in (|4.1| ) have the same scale dependence. Thus, we 
recover the standard prediction for the scalar spectral tilt in single field inflation. The 
function /2(C) which appears in ( |4.1| ) is plotted in Fig. 0. For £ > 1 eq. (|1|) provides 
an asymptotic expression for f'2 while, on the other hand, eq. (|3.20| ) provides a good fit 
in the 2 < £ < 3 interval (as we discuss below, this is the most interesting interval for 
phenomenology). The COBE normalization PAk) = 25 • 10 -10 is satisfied along the curve 



^COBE (£) 



-4tt£ 



l + A /l + 10-8/ 2 (Oe 4 ^ 



(4.3) 



2/2(0 

in the £ — V plane. We see that, at low £, the contribution to ( [4.1| ) sourced by the gauge 
field is subdominant and we recover the standard result Vqq-q^ — 5 • 10 -5 . The two 

contributions in ( [4.1[ ) become equal at £ = 2.9. As £ is increased, inverse decay effects 

— 1/2 
dominate the spectrum and the value of "P^obe mus t be (exponentially) decreased to avoid 

over-producing density perturbations. 
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Figure 4: Values of parameters leading to the observed COBE normalization of the power spectrum 
(red line), and reference values for the nongaussianity parameter f^* 1 = 10,266,8000 along this 
curve. See the main text for details. 



The curve ( [4.3[ ) is shown (red solid line) in Figure The black dashed line shown in 
the Figure separates the (lower) region in which the vacuum fluctuations dominate from 
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the (upper) region in which the fluctuations from inverse decay dominate. In the region 
above the solid black line in Fig. |] the backreaction bound (2.16) is violated. In that 



region of parameter space, the production of gauge field fluctuations is so strong that 
dissipative effects (rather than the potential V ((/))) dominate the motion of 0. From Fig. [I] 
we see that this backreaction effect can be safely disregarded after the COBE normalization 
is imposed. 7 Finally, Figure |I| shows some reference values of /jy^ 11 along the COBE 



normalized curve; we discuss this in subsection 4.2. 
4.2 Tensor-to-Scalar Ratio 



Similarly to the result fl4.1| ), the GW spectrum is also the sum of two uncorrelated contri- 
butions, one due the modes generated from the vacuum, and one due to the modes sources 
by the gauge field quanta. In subsection l3~5| we derived the result 



2 / u \ n T 



Pgw 



2H 2 f k 



H 



2 



1 + ^(0 



(4.4) 



7r 2 M p 2 \k 

where tit = — 2e. The second term in the square braces corresponds to the gravitational 
waves sourced by gauge field quanta. The function f^L is plotted in Figure 0. Eq. (3.41) 



provides a large argument expansion while eq. fl3-42| ) provides a fit in the 2 < £ < 3 interval 
(as we discuss below, this is the most interesting interval for phenomenology). 

We define the tensor-to-scalar ratio in the usual way, by normalizing the amplitude of 
the power in gravitational waves to that in scalar fluctuations 

r = —— (4.5) 

For £ < 1, inverse decay effects are negligible and we recover the standard consistency 
relation for r, familiar from single field inflation. At £ — > oo, on the other hand, r tends 
to a different constant value, which is smaller than the asyptotic £ — >• value. From (tO), 
Q4.4j), and from the result tit = — 2e it is easy to show that 



-8n T if £ < 3; 
1.8 ny as £ — > oo. 



2 . (4-6) 



Therefore, r is independent of £ if either vacuum fluctuations or inverse decay effects 
dominate, while it interpolates between these two asymptotic values for intermediate £. 

4.3 Nongaussianity 

As discussed in section ||, the cosmological fluctuations generated by inverse decay effects 
are highly nongaussian. There are many different ways to parametrize departures from 
gaussianity. A standard work- horse is the local ansatz: 

C(x) = ( g (x) + ^/jvT 1 [e g (x) - (C 2 g (x))] (4.7) 



7 One should also ensure that the energy density of the produced quanta gives a negligible contribution 
to the expansion of the universe. The resulting condition, eq. ( |2.17| ), cannot be shown in the £ — V plane, 
and therefore needs to be studied on a case by case basis. We have verified that this condition is satisfied 
for the models studied in this section. 
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where Q g is a gaussian random field and f^f^ quantifies the amount of nongaussianity. 
Although this simple parametrization has received considerable attention, it is certainly not 
the only well-motivated model for a nongaussian curvature perturbation. More generally, 
one should consider the bispectrum, B^(ki), which is the 3-point correlation function of £ 
in Fourier space: 

<C kl Ck 2 Ck3> = B ( (ki) 6® (k x + k 2 + k 3 ) (4.8) 

The bispectrum is a function of three momenta, ki, that form a triangle: ki + k 2 + 
k 3 = 0. Hence a generic bispectrum may be characterized by specifying three interesting 
properties: the magnitude of the function, its dependence on the shape of the triangle, 
and its dependence of the size of the triangle. These properties are usually referred to 
as the size, shape and running of the nongaussianity, respectively. To characterize these 
properties, we find it convenient to introduce 

|ki| = A; , \h 2 \=x 2 k , |k 3 | = X3 A; (4.9) 

so that k encodes the overall size of the triangle while the dimensionless quantities x 2 , £3 
encode its shape. 

If we assume a local ansatz ( |4.7| ) then the bispectrum has a very particular dependence 
on momenta: 



?J ocal (^) = ^(2nf 2 f%fP ( (k) 2 ^ 



3(2vr) 5 / 2 P C (fc) 2 1+4 + 4 f local u , m 

- ~Io ¥ 4xY~ tNL { ' 

This function peaks in the squeezed limit where one of the wave-numbers is much smaller 
than the other two (for example k\ <C k 2 ,ks). 

4.3.1 The Size and Running of the Nongaussianity 

The bispectrum from axion inflation contains two uncorrelated contributions correspond- 
ing, respectively, to the usual vacuum fluctuations and to the fluctuations generated by 
inverse decay processes. As is well known, the former contribution gives rise to unde- 
tectably small nongaussianity and may be ignored. The second contribution, however, is 
more interesting and this was computed in subsection [OJ. We found that the bispectrum 
from axion inflation is very different from the local form ( [4. 10 ). The bispectrum from 



axion inflation peaks on equilateral, rather than squeezed, triangles. Nevertheless, it is 
conventional to characterize the size of nongaussianity by matching to ( [4.10 ) on equilateral 
triangles |ki| = | ki2 1 = | IC3 1 . Proceeding in this way we find 

fequi i_ Mjj h *) V ^ G< am 
/nl - pdk) 2 (4- 11 ) 

The function / 3 (£; 1, 1) is plotted in Figure Eq. (|3.32|) provides a large argument 
expansion of this function, while eq. ( |3.33D provides a fit in the 2 < £ < 3 interval. 
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Figure 5: Value of (£) due to the perturbations £ sourced by the gauge quanta. Also shown 
is the 95% CL upper bound from WMAP 7 [M. 



We stress that this result does not include the negligible contribution from 
accurate as long as l/^ 11 ) 1. 



and is 



From eqs. Q4.3| ) and ( 4. lip , from the observed power P^ = 25 • 10 -10 (we disregard 
the slow-roll suppressed scale dependence of Jnl), and from the expressions of f 2 and 
computed in the previous section, we can compute the value of f^i (0 along the COBE 
normalized curve. We show this in Figure || (see also the reference values shown in Figure 
H|). We notice that f^ 11 (£) saturates to a constant value at large £ (in the region where 
the vacuum contribution to P^ is negligible). From the large value asymptotic expressions 
(|3l9|) and (|^3|), we find (oo) ^ 8,600. This value is already above the 95% CL 



upper -214 < /° q f uil < 266 obtained from the WMAP 7 data |39(]. This limit rules out 



> 



2.65. The Planck satellite is expected to be able to resolve to O(10). 

Having quantified the size of nongaussianity in axion inflation, we now turn our atten- 
tion to its running. From eq. ( [3.28 ) we can see that 



B i (k 1 ,k 2 ,k 3 ) = k B i (l,x 2 ,x 3 



(4.12) 



disregarding the mild, slow roll suppressed scale dependence of the vacuum solutions. The 
overall A;" 6 behavior reflects the near scale invariance of the bispectrum from axion inflation. 



o!ln|/jVLl 
dink 



Slight departures from scale invariance are quantified by the index ung — 1 
which is easily seen to be proportional to the slow roll parameters e, 77, and hence negligible 
whenever the observational bounds on n s are satisfied. We conclude that the running of 
nongaussianity is uninterestingly small in axion inflation. 



4.3.2 The Shape of the Nongaussianity 

In order to discuss the shape of the bispectrum, it is natural to extract the strong A;~ 6 
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scaling in ( |4.12| ) and define a "shape function" of the form 

S(ki) = NihhhfB^h) 



(4.13) 



where the constant of proportionality, N, is arbitrary. This shape function coincides with 
the quantity that was plotted in many previous works, including |47]] for example. For the 
case of interest, we have 



S (€;x 2 , x 3 ) 



l + xl + xl / 3 (£; x 2 , x 3 ) 



X2X3 3/3 (£; I, 1) 



(4.14) 



which is normalized so that 5(1, 1) = 1. Note that the bispectrum is defined only in the 
region x 2 + x% > 1, which follows from the triangle inequality. Moreover, the bispectrum 
is symmetric under interchange of any two momenta, and therefore we can restrict to the 
region X3 < x 2 < 1 to avoid considering the same configuration more than once. 

We plot the shape function S(x 2 , £3) from axion inflation in the left panel of Fig. |6|. 
The bispectrum in this model depends on the parameter £. In practice, however, we find 
that only the size of the nongaussianity (quantified by ff^ 11 ) depends strongly on £. The 
shape function S(x 2 ,x 3 ), on the other hand, is very mildly dependent on £. In Fig. || 
we work in the £ — > 00 limit, in which case the shape becomes independent of model 
parameters. (This can be seen by using the large argument expansion ( |3.18| ) of X in the 
expression for / 3 .) For £ ~ 0(1) this fi gure would be nearly indistinguishable. 
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Figure 6: In the left panel we plot the shape function S (X2, £3) in axion inflation, showing that 
this peaks on equilateral triangles. We work in the limit £ — ¥ 00, however, this Figure would be 
nearly indistinguishable had we chosen £ = 0(1). In the right panel, for comparison, we plot the 
analogous shape function obtained from the standard equilateral template. 

From Fig. || we see that the bispectrum from axion inflation peaks on equilateral 
triangles (corresponding to x 2 = X3 = 1) and is thus qualitatively similar to the so-called 
equilateral template which is often employed to analyze CMB data [48, ^] 



B, 



cqm\{h) OC ,3,3 



1.31,3 
1"^3 



1.31.3 
K 2 ^3 



h 2 h 2 h 2 



+ 



k\ 



+ (5 perms) 



(4.15) 
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(where the permutations only act on the last term). Equation ( 4.15| ) is the template that 



is used to obtain the WMAP7 limit on /nl" 1 ^ ^66 which we employed in subsection |4.3.1 
The shape function associated with the template ( 4.15| ) is plotted in the right panel of 
Fig. U for comparison with the analogous result from axion inflation. 

To quantitatively compare the bispectrum from axion inflation to the equilateral tem- 
plate ( 4.15| ) we follow [47 1 and define a scalar product between any two bi-spectra as 

Bi ■ B 2 = ^ B 1 (k 1 ,k 2 ,k 3 )B 2 (k 1 ,k 2 ,k 3 )/ ««) (4.16) 
where o\. is the variance of a given mode and the summation runs over all possible trian- 



gles. As shown in [47], this product is the best estimator for the overlapping of any two 
distributions: if we assume that the real data have the bi-spectrum Bi, the template B 2 
will produce a higher / lower value of nongaussianity according to how large / small the 
product fl4.16|) is 47]. To be quantitative, one defines the cosine of the "angle" between 



the two bi-spectra as []47| 



cos (Si, B 2 ) = 



B\ ■ B 2 



(B 1 .B 1 y/^B 2 -B 2 y/^ 



(4.17) 



The cosine can be used to quickly estimate how well the limit given in the literature on 
some given template applies to a different shape [47]. 

We have computed the cosine of the "angle" between the bispectrum from axion infla- 
tion and the equilateral template for several different values of £. These results are reported 
in Table 1. There we see that the cosine depends only very weakly on £ and saturates to a 
value = 0.93 in the limit £ —> oo. This confirms quantitatively our previous claim that the 
shape of nongaussianity is insensitive to £. 





COS (-Bjnv. decay ) -^equil) 


COS (-Binv. decay j -Berth) 


2 


0.94 


-0.093 


3 


0.94 


-0.12 


5 


0.93 


-0.13 


oo 


0.93 


-0.15 



Table 1: Cosine of the "Overlapping angle", eq. (1.17), between the nongaussian shape generated by the inverse 



decay, and the equilateral (column 2) and orthogonal (column 3) templates, for different values of £. 



Table 1 shows that cos (-Bmv.decayj B equ ii) is very close to unity and hence we expect 
that the WMAP7 limit -214 < /^f < 266 can be applied also to axion inflation (at 
least to first approximation). This justifies our interpretation of the observational limit on 
nongaussianity in subsection [4.3.1 . 

Although the bispectrum from axion inflation is very similar to the equilateral tem- 
plate, the two shapes are not identical. It may be interesting to characterize the difference 
between these two shapes - indeed, this would become pressing in the event that Planck, or 
some other future mission, should detect a non-vanishing bispectrum on equilateral trian- 
gles. From Figure ^, we see that the bispectrum from inverse decay mostly differs from the 
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equilateral template for X2 — £3 — 1/2, corresponding to "flattened" triangles, where one 
side is twice the length of the other two. For such triangles, the bispectrum from inverse 
decay is significantly greater than the equilateral template. Hence, this provides a natural 
limit in which the two shapes may be distinguished. 

The fact that axion inflation gives a large nongaussianity on flattened triangles implies 
a small but nontrivial overlap with the so-called "orthogonal template" that was introduced 
in @: 

3 3 3 8 3 

B rth(h) ^-T^-pU~T3p- T2T2T2 + T7E2T3 + ( 5 P ermS ) ( 4 - 18 ) 
1 2 1 3 2 3 1 2 3 1 9 3 

(where the permutations only act on the last term). The corresponding shape function 
evaluates to +1 in the equilateral limit, and to —2 along the X3 = 1 — X2 boundary 
(which includes the flatten triangle configurations). In Table 1 we have computed the 



cosine of the angle between the bispectrum from axion inflation and the template ( 4.18| ). 
We find cos (i?i nv .decay 5 -Borth) — —0.15 (at large £, see Table 1 for intermediate £), while 
cos (5 equi i, 5 ort h) = 0.21. Therefore, the overlap with the orthogonal template may pro- 
vide a useful tool to discriminate observationally between the nongaussianity from axion 
inflation and that of the equilateral template. 

4.4 Large Field Inflation 

Once the COBE normalization is imposed, the key phenomenological predictions of any 
inflationary model are the spectral index, n s , the tensor-to-scalar ratio r, and the nonlin- 
earity parameter Jnl- As we discussed previously, the spectral index in axion inflation 
has the standard form n s = 1 + 2r/ — 6e and requires no further discussion. The remaining 
observables depend on the coupling a/f, the axion velocity (p and the Hubble rate H. Out 
of these three quantities, we have defined the two combinations 

e = p i/2 = H_ u 19) 

5 ~ 2fH ' ' ~ 2w\<p\ { ' 

Both the two and three point correlation functions of £ can be written solely in terms of 
these two combinations. Therefore both the power spectrum and the bi-spectrum are a 
functions of £ and V only. The COBE normalization fixes V 1 / 2 in terms of £, see Figure 
|I[ Therefore, the predicted nonlinearity parameter /jvl is function of £ only, see Figure |5|. 

The tensor to scalar ratio, on the other hand, is a function of a different combination 
of parameters, see eq. ( |4,4j ). For this reason, we cannot present it in a plot as a function 
of £ only. This can only be achieved once the potential for the axion is specified, since this 
provides one additional (slow-roll) relation between the parameters of the model. 

For a specific choice of inflationary potential V((f>), it is of interest to determine how 
the combinations ( 4.19p , along with the tensor-to-scalar ratio, depend on parameters of 



the underlying theory. As we discuss in [|T^] and also in Section ^ in the models of 
axion inflation of interest, the axion/inflaton dynamics effectively occurs as in a large field 
inflationary model with potential: 

V = —(fP (4.20) 
P 
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where X„ has mass dimension 4 — p. Using slow roll approximation, we find 



\3p M i-p/2 



— hm v , e 

2N p ' s 



^7" 



-4- (4.21) 



where Af is the number of e-folds between the moment the CMB scales left the horizon and 
the end of inflation. 
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Figure 7: Predicted values for the equilateral Jnl parameter, and for the tensor to scalar ratio r 
in axion inflation models, as a function of the coupling of the axion to gauge fields, when the axion 
dynamics is effectively described by a potential V oc ip p , with p = 1,2. 

The value of \ p is fixed by COBE normalization and we assume N = 60 e-foldings 
of inflation. Then, for any given value of p, all observational predictions can be written 
in terms of //a only. In axion monodromy [12], p = 1; in most of the other models one 
expands the potential close to the minimum, where it is quadratic, p = 2. We therefore 
show in Figure |7] the predicted values of Jnl and of r for these two cases. We notice that, 
once /jvl is required to be below the WMAP7 bound, the standard value for r is recovered. 

It is interesting to note that axion inflation models generically predict the same values 
of n s and r as would be obtained in vanilla chaotic inflation. However, our scenario predicts 
also a large nongaussianity with a (nearly) equilateral shape. Axion inflation provides a 
rare example of a theory which predicts both a detectable tensor-to-scalar ratio and a large 
equilateral bispectrum. Note that, if such a nongaussian signal is eventually detected, then 
it will immediately fix the value of the coupling a/f. On the other hand, if Planck fails 
to detect nongaussianity then we will have a surprisingly stringent bound on the strongest 
axion-type couplings between the inflaton and any gauge field. 

5. Cosmological Perturbation Theory 

In section ^ we summarized the results of ]l8[, providing a brief overview of the growth 
of gauge quanta and the production of inflaton fluctuations via inverse decay processes, 
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in the model (2.1). In section |3] we provided a detailed computation of the relevant cor- 



relation functions. There, and also in [18|, we neglected scalar metric perturbations. It 
is intuitively clear that this should provide a sensible approximation when /jvx ~ 1> since 
the nongaussianity associated with gravitational interactions is tiny as compared to the 
axion-type interactions that we study (1.3). In this section, we reconsider the analysis of 
section |3|, this time consistently accounting for scalar metric perturbations. As would be 
expected, there will be no significant change in our key results. 

5.1 ADM Formalism 



We consider the theory (2.1), minimally coupled to Einstein gravity 

1„ . 1 a 



S 



d*Xy 



-g 



Ml 



1 [iv 1 ap 



(5.1) 



The Levi-Civita tensor e^ va ^ is defined as 



-V 



-g 



(5.2) 



where the alternating symbol Tf- va ^ is +1 for even permutations of its indices, —1 for odd 
permutations, and zero otherwise. 

To study the cosmological perturbations it is convenient to employ the Arnowitt-Deser- 
Misner (ADM) form of the metric 



ds 2 



N 2 dt 2 + hij(dx l + N l dt)(dx 3 + N'dt) 



(5.3) 



This parametrization has the advantage that the lapse function, N, and shift vector, N l , 
appear in the action as Lagrange multipliers and hence are not dynamical degrees of free- 
dom. We consider only scalar perturbations of the metric, and we employ the following 
gauge-fixing choices 



<p(t,x) = 4>{t) + 5<p(t,x) 
hij(t,yL) = a 2 (t)6ij 
A = 



(5.4) 
(5.5) 
(5.6) 



corresponding to Coulomb gauge for the vector, along with a flat slicing of the space- 
time (since the linearized solutions are used for A^, we can also impose diAi = 0). It is 
now straightforward to vary the action ( |5.1[ ) with respect to N, N l , yielding the constraint 
equations. These constraints may be solved perturbatively so that N, N l can be eliminated 
in favor of 5(p, Ai. In this manner, we can derive an action which describes the dynamics of 
the scalar and gauge field fluctuations, and their leading interactions. This procedure has 



been described many times in the literature, following [21|, and we do not report the explicit 
steps of the computation in full detail in this work. We also note that the cosmological 
perturbations in a model very similar to ours (differing only by the inclusion of the term 
/ (ip) F 2 rather than ipFF) have already been studied in detail in [40|. 
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5.2 Quadratic Action 



The first non-trivial contribution in the perturbative expansion of ( p,l[) is the quadratic 
action for the fluctuations dtp, A$, which give rise to the linearized equations of motion. To 
leading order in slow roll parameters, we find the following result for the quadratic terms 
involving the scalar field fluctuation 



St 



hi 



drd xa 



(d T 5tp) 2 - diS^didv + (a 2 m 2 - 3-^ j 6<p 2 



(5.7) 



The dynamics of the gauge field, on the other hand, arise from the quadratic terms in the 
expansion of the Maxwell term in the action Q5.1|) . We find 



- / drd 3 x 
2 



A>A> - d,A ldj A t + ^e ijk A i9j A k ^ 



(5.8) 



5.3 Cubic Action 

The leading nongaussian effects are encoded by the cubic terms in the expansion of (|5.1|), 
These may be divided into two categories: interactions involving 3 scalars, and those 
involving one scalar and 2 vectors. In the first category, we find 



s 3 



/ 



,2 jJ 



drd 6 x 



W.M 2 



[26ip'did- 2 (6<p')di(6ip) - MV) 2 - <fy>$ 0^0] ( 5 - 9 ) 



where we have introduced the inverse Laplacian operator, defined through the relation 
d~ 2 di&if = f. As expected, this coincides with the well-known result for single field 
inflation [fTj], ||. In Q we have neglected a term proportional to V'" , which is usually 
sub-dominant in the slow roll expansion (see [^] for more discussion). 

The interactions ([^) would be present even in the absence of the gauge field A^. The 
contribution to the non-linearity parameter due to these terms is proportional to slow-roll 
parameters pl|, 22, E3| and may therefore be neglected whenever f^L ~ C(l). In this 



work, we are therefore justified to ignore (p.9|) . 

Much more interesting for our analysis is the second category of interaction: those 
involving one scalar 5<p and two gauge fields Ai. From the expansion of the Maxwell term 
F 2 in the action (|5.1|) we find 



Sf A D 



drd s x 



S(p 



(5.10) 



which agrees with [40]. Finally, from the last term in ( |5.1[ ) - the pseudo-scalar interaction 
- we have a contribution 



St AA D J drd 3 x j 5p [-e^A'&A, 



(5.11) 



This interaction gives rise to the source term in ( 2.18| ) and was accounted for by the analysis 
of [ 1 8 1 . Equations (5.9, 5.10| , 5.11| ) exhaust all interactions at cubic order. 
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The interaction ( 5.11 ) would be present even in the absence of metric fluctuations and 
was already accounted for in the analysis of fl~8[| and section ||; this term simply arises from 
J1.3D . On the other hand, the interactions ( |5.9| ) and ( 5.10| ) are "new" and arise due to 
the consistent inclusion of metric perturbations. It is easy to see already that these "new" 
interactions are completely negligible as compared to the pseudo-scalar coupling ( 5.11| ). To 
see this, note that the strength of the interaction ( |5.11 ) is controlled by the axion decay 
constant which is 

- > 10 2 — 

f~ M, 

whenever the nongaussianity is observationally interesting [|j 
strength associated with the interaction ( 5.10| ) is instead 



(5.12) 

The analogous coupling 



10 



-i 



1 



(5.13) 



where in the last approximation we have assumed n s — 1 ~ e ~ 10~ 2 , which is true in most 
axion inflation models. A similar line of reasoning can be used to verify that the trilinear 
interactions described by ( |5.9| ) are slow roll suppressed with respect to M~ l and are known 
to yield a negligible contribution to Jnl [ 43 1 , as compared to (|5.1l| ). 

This analysis indicates that the inclusion of metric perturbations has no significant im- 
pact in the interesting regime where Jnl ~ 0(1)- This intuitive idea has been exploited on 
numerous occasions in the literature and was recently conjectured as a general decoupling 



principle in (44]. 



5.4 The Equations of Motion 



In subsections p.2| and |5.3| we derived the leading order contributions to the action for 
the scalar and gauge field fluctuations. In order to compute the nongaussianity in the 
model, there are two main approaches. We could use the in-in formalism to compute the 
bispectrum directly from the interactions ( f>$ , [5,10 , 5.11| ). This is the approach that was 



adopted in [21, 22, 23, 40, 41]. Alternatively, we can use the perturbed action to derive the 
Klein-Gordon equation, and then compute the bispectrum following ^3|. Both approaches 
give the same answer. In our case, we find the second approach to be more convenient, 



since it will allow us to take maximal advantage of our previous analysis [18]. 

Variation of the quadratic action ( |5.8| ) gives the linear equation of motion of the gauge 

field 



o < 



8;- V - y VX 



A(t,x) 







(5.14) 



This equation coincides exactly with our previous result ( |2.10| ). The consistent inclusion 
of metric perturbations has no impact on the linear fluctuations of the gauge field, which 
is a consequence of the fact that (Aj) = 0. 

Next, we consider the Klein-Gordon equation for the scalar fluctuation 5<p. Variation 
of ( |5.7| ) along with (|5T9,5.1C,5.11) gives rise to the following equation of motion 

..aV 

d 2 + 2Ud T 



V 2 + 



2 2 

a m 



3<P 



M 2 
p 



o 



6<p = ja 2 E ■ B 
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,2i/ 



2M| 



p2 , r>2 i 

t + ^a- 2 9 r ( a 4 V • ( £ x B 



+ 



(5.15) 



where the • • • on the last line denotes terms of order (5ip) 2 which arise due to variation of 
(15.91) . As discussed in subsection 15. 3L those terms contribute negligibly to the non-linearity 



parameter and can be ignored. As a check on our results, we have also derived equation 



(|5.15|) using a completely independent analysis: by expanding the Klein-Gordon equation 
to second order in perturbation theory and using the Einstein constraint equations to 
eliminate the metric fluctuations, following [45|. This analysis is briefly summarized in 
Appendix B. 

Equation ( |5.15|) differs from the result ( |2.18| ), that was derived by neglecting metric 
perturbations, in two respects. First off, equation ( |5.15| ) includes a slow-roll suppressed 
correction to the effective mass of the inflaton. The second discrepancy is the additional 
contribution to the source on the second line of ( 5.15| ), proportional to M~ 2 , arising from 
variation of ( |5.10 ). This contribution is parametrically of order ~ V^ E m B an d may be 



neglected as compared to the term jE ■ B, as discussed in subsection (|5.3|) . The results 
of this subsection provide rigorous justification for our previous use of equation ( p. 18 ) to 
study the cosmological fluctuations in the model (|5.1[) , in the /tvl ~ 1 regime. 



5.5 The Curvature Perturbation 

It now remains to construct the curvature fluctuation on uniform density hyper surf aces, 
C- Beyond linear order, the relationship between £ and Sip is simplest to derive using the 



S N- formalism |46|]. Because the gauge field does not contribute to the expansion history 
of the universe, we have the simple result on large scales: 



C = 5N 



ON . Id 2 AT , 2 



+ ... 



-—Sip 

<P 

-—Sip 



l + (2e 



(5.16) 



where on the last line we have neglected a slow-roll suppressed correction which contributes 
negligibly to Jnl in the case of interest. 

Equations fl5.15| ) and ( 5.1 6| ) are the main results of this section. Taken together, they 



provide a rigorous justification for the analysis of IS], and of the previous sections, which 
disregard the effect of metric perturbations. 



6. Models of Axion Inflation 



Our results concerning nongaussianity and inverse decay effects are quite general and may 
be applied directly to a variety axion inflationary models. Moreover, we expect that our 
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qualitative results may have also implications more broadly, for example in any multi- 
field scenario that involves a dynamical axion. 8 In this section we survey some interesting 
microscopic constructions, both from field theory and string theory, and comment on the 
possibility of large nongaussianity. The key input which our analysis requires from a 
microscopic computation are the inflaton potential V(ip), the decay constant /, and the 
dimensionless parameter a. Together, these determine the quantity 

aMl V' 



2/ 



V 



(6.1) 



that measures the strength of inverse decay effects. As we have seen, inverse decay processes 
cannot be neglected when £ > 0(1) which roughly translated into f /(aM p ) <, 10 -2 for the 
most interesting models; see Figure [7|. Whenever this inequality is satisfied, our findings 
strongly affect the phenomenology of these models. 

In most of the scenarios that we survey, both V(<p) and / tend to be fairly well 
understood while, on the other hand, the coefficient a is rather more model dependent 
(although calculable in principle). A detailed computation of a in each interesting scenario 
is beyond the scope of this paper and we leave such an analysis to future work. Although we 
generically expect a = 0(1), it should nevertheless be noted that this parameter contributes 
a source of theoretical uncertainty to what follows. We believe that our phenomenological 
results should provide a motivation for a detailed microscopic computation of a in the 
various scenarios discussed below. 

6.1 Natural Inflation 

The original natural inflation proposal || was based on the potential 9 



V(<p) A 4 



1 — cos 



(6.2) 



For a spectral index n s ^ 0.95 this model requires / £ \^8tt M p [16], for which the inflaton 
dynamics shows little difference from that standard chaotic inflation with V(ip) = ^m 2 ip 2 , 
at least for the present considerations. Such large values of / weaken the pseudo-scalar 
coupling cpFF and inverse decay effects are negligible unless a £ few x 10 2 . While specific 
situations have been constructed that can result in a large a - for instance in the some 



extra-dimensional model [15] - such large values conflicts with our general expectation that 
a = 0(1). We conclude that large nongaussianity seems unlikely in the simplest models of 
natural inflation. 

In spite of apparent simplicity, however, the original natural inflation model H seems 
incompatible with UV completion. If we interpret <p as a PNGB then / > M p suggests 



8 This seems to be a generic expectation for closed string inflation models. In the effective SUGRA 
description the Kahler moduli n are typically paired with axions 0i into complex fields Ti = Ti + iQi. Absent 
tuning, one expects the curvature of the scalar potential to be comparable in the n and 8i directions. Hence, 
nontrivial modular dynamics is typically also associated with nontrivial dynamics in the axion sector. See, 
for example, reference ^ij] for a review. 

9 See for some proposed particle physics realizations. In E3] this model was realized in supergravity 
- the so-called "axion valley" model. 
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a global symmetry broken above the quantum gravity scale, where effective field theory 
is presumably not valid. Moreover, / > M p does not seem possible in a controlled limit 
of string theory [fj"7f | (which is the only known framework wherein such questions may 
be addressed). Hence, requiring the existence of a sensible UV completion automatically 
pushes us towards the regime / <C M p where inverse decay effects are important. We will 
illustrate how this works in several explicit microscopic realizations below, however, it is 
clear that this trend applies more generally. 

6.2 Double-Axion Inflation 

Perhaps the simplest scenario to realize natural inflation with / < M p is the double-axion 
model proposed in ||. This model is characterized by two axions, 6 and p, whose potential 



1 - cos ( — + — 

h 9i 



(6.3) 



i=i 

arises from pseudo-scalar couplings to two different gauge groups: 4-FiFi and —FiFi. For 

Ji 9i 

f\92 = 51/2 one linear combination of 9 and p becomes a flat direction of the potential, 
corresponding to an enhanced symmetry of the theory. Taking f\g2 — 31/2 the curvature in 
this direction becomes controllably flat and one may realize natural inflation with fi,gi <C 
M p . This scenario therefore provides a simple and compelling illustration of how the 
requirement of a UV completion leads to large nongaussianity in natural inflation. 

6.3 N-flation 

The N-flation model |i~0[ ] is based on N axion fields (pi, each with its own softly broken shift 
symmetry, resulting in a separable potential of the form 

= E A * I 1 - cos (f ) ] = E ( 6 - 4 ) 

where, in the second equality, we have expanded in small field values (pi -C /j. For JV S> 1 



the assisted inflation mechanism 51] allows inflation to proceed even while all the decay 
constants are sub-Planckian, /j < M p . To a first approximation, the dynamics may be cap- 
tured by a simple quadratic potential V e g = ^m 2 <3? 2 for the collective field 3> 2 = £^ 93? 1 11 



Successful inflation requires that the collective field traverses a super-Planckian distance 
in field space, A$ ■> M p . This is achievable with sub-Planckian ipi provided the number 
of axions is sufficiently large. Roughly speaking we require p 



N ~ 240 ( -j^- ) (6.5) 



/ avg 



For / avg ~ 10 _1 M P we have N ~ 10 3 while / avg ~ 10 _2 M p would require N ~ 10 6 . 

N-flation makes sense as a purely field theoretical construction, however, much of the 
interest in this scenario arises because the exponentially large values of N that are required 
may be rather generic in string theory compactifications. In this context, axions arise on 
dimensional reduction, from integrating p-form gauge potentials over p-cycles (see [4S, 53] 



- 32 - 



for a review). For example, in type IIB string theory one has axions hi and Cj which 
arise, respectively, from integrating the Neveu-Schwarz (NS) and Ramond-Ramond (RR) 
2-forms Bmn an d Cmn over compact 2-cycles. Generically, instanton effects break the 
shift symmetry down to a subgroup 6j — > bi + (27r)/j, leading to periodic contributions to 
the effective potential (similarly for c$). 

The general framework described above shows how N-flation may arise within string 
theory: the low energy theory contains one axion for each independent cycle that Bmn can 
wrap and generic Calabi-Yau compactifications may contain exponentially large numbers 
of such cycles. There exist know examples with N as large as ~ 10 5 [plj , however, to our 
knowledge there is no general theorem that prohibits finding Calabi-Yau manifolds with 
even larger values. 

There may be many ways to realize N-flation in an explicit, stabilized string theory 



compactification. The first efforts in this direction were undertaken in [11 1 and further 
analyzed in [52]. See also [55] for an alternative construction. 



In order to quantify inverse decay effects, we are most interested in the effective cou- 
pling between the collective field 3> and a given gauge field. The detailed derivation of such 



interactions is rather complicated and dependent on model building details (see [56, 57, 58 
for a more detailed discussion). We will not attempt to estimate the effective value of 
a//, but rather note that a generic spectrum is expected to contain some axions with 
fi <C M p , hence strong inverse decay effects (and their associated nongaussianities) are at 
least plausible in N-flation. 

6.4 Axion Monodromy Inflation 

The axion monodromy model [12| is a string theoretic construction based on a single 



axion field. The key ingredient in this construction is a suitably wrapped brane which 
leads to a non-periodic contribution to the inflaton potential, explicitly breaking the shift 
symmetry. This monodromy in the moduli space allows the axion to develop a kinematically 
unbounded field range and accommodate the super-Planckian excursions required for large- 
field inflation. In the effective field theory description of this scenario, the inflaton potential 
has the form 

V(Lp)=» 3 <p + A 4 cos(^) (6.6) 

where the linear term arises from wrapped branes while the (subdominant) periodic mod- 
ulation arises from instanton effects. 



In [13] the decay constant for axion monodromy inflation was studied in detail. It was 



shown that microphysics bounds the allowed values as 

1/4 , 

°- 06 ^ < 4 <0 - 9 ^ (6 - 7) 

where g s < 1 is the string coupling and V S> 1 is the compactification volume in string 
units. This bound illustrates that / <C M p in this model. If a = 0(1) can be realized in 
a consistent string compactification, then large nongaussianity is easily accommodated by 
axion monodromy inflation. 
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The small periodic modulation of the potential ( p,6|) is not important for determining 
the number of e-foldings of inflation, however, this term may nevertheless have an important 
role of the phenomenology of the model. In 12], and also in [59], it was shown that axion 



monodromy models can give rise to large resonant-type nongaussianities [30], 44]. 

Depending on the parameters, either resonant effects or inverse decay effects may dominate 
the bispectrum in the regime / <C M p . It would be interesting to study the combined 
observational impact of these two effects in future work. 

6.5 Axion/4-Form Mixing 

In [|14|] a realization of chaotic inflation was proposed which shares many features of the 
axion monodromy model discussed above. Here the axion (p mixes with a 4-form field 
strength through a term like ip e ,J,ua ^ F^uaP ■ The theory also includes charged membranes 
which source a background for the conjugate momenta of the gauge field and break the 
axion shift symmetry, giving rise to a non-periodic potential V(ip) = ^m 2 ip 2 which is robust 



against a wide variety of corrections [61]. There is no obstruction to realizing inflation with 
/ <C M p in this model, and hence we generically expect that inverse decay processes may 
play an important role. 

As in the case of axion monodromy, the power-law potential will generically be modu- 
lated by subdominant oscillatory features arising from instanton effects. In IHJ nongaus- 
sianities were discussed. In general resonant effects may operate in concert with inverse 
decay processes, and it would be interesting to study their combined impact. 

6.6 Dante's Inferno 



In |62j a model was proposed which consists of two axions, r and 9, with decay constants 
fr < fe M p . It is assumed that a linear combination of these receives a periodic potential 
from nonperturbative effects and, moreover, that some explicit shift symmetry breaking 
effect generates a non-periodic contribution W[r). The effective potential therefore takes 
the form 

/ r 

. fr fe 



V(r,9) = W(r) + A 4 



1 — cos 



(6.8) 



For a power-law W(r) = /i 4_p r p the dynamics of this model are well approximated by a 
single field with effective potential V e R = {fr/ fe) p ^~ p f P c s anc ^ our analysis of inverse decay 
processes is directly applicable. The possibility of large nongaussianity follows immediately 
from f r < / e < M p . 



In |52] it was shown how to embed the model (6.8) within string theory as a modest 



extension of axion monodromy [12, 13 1. It may be possible to realize such a potential from 



axion/4-form mixing, along the lines of [14, 61]. 
6.7 Multi-Field Scenarios 

So far we have focused our attention on models which can, at least to first approximation, 
be described in terms of the dynamics of a single field (that may represent a collective 
excitation or be related non-trivially to the axions of the original theory). Such scenarios 
are appealing, because our analysis of inverse decay processes applies more-or-less without 
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modification. However, there are also a number of interesting inflation models which in- 
volve dynamical axions but are not well described in terms of a single dynamical field. (As 
discussed in footnote ||, this scenario seems especially natural in supergravity models.) In 
multi-field models with / <C M p strong inverse decay effects are clearly possible, however, 
the detailed phenomenology is more complicated than what we have presented. Neverthe- 
less, we expect that large nongaussianity should be possible. It would be interesting to 
explore multi-field inverse decay processes in future works. 

There are a variety of interesting multi-field models involving axions, for example 



racetrack inflation [63, |6_|] and the axionic D3/D7 model [65]. Possibly, interesting models 



from the perspective of obtaining strong inverse decay effects, are based on the large volume 



compactification of |66|, |67], j^g]. For example, roulette inflation [S9| — which generalizes [70| 



to incorporate the dynamics of the axionic partner of the Kahler modulus - is characterized 



by significant motion in the axion direction during the early stages of inflation. See [71, 72] 
for a discussion of axion/moduli couplings and decays in large volume inflationary scenarios. 
See also [^] for a related model. 

7. Conclusions 

Probably the greatest difficulty in inflationary model-building is to protect the required 
flatness of the potential from radiative corrections. One of the simplest solutions to this 
problem is to assume that the inflaton ip is a Pseudo-Nambu-Goldstone-Boson (PNGB). 
In this case the inflaton enjoys a shift symmetry ip — > ip + const, which is broken either 
explicitly or by quantum effects. In the limit of exact symmetry, the ip direction is flat 
and thus loop corrections to the potential are controlled by the smallness of the symmetry 
breaking. PNGBs, like the axion, are ubiquitous in particle physics: they arise whenever 
an approximate global symmetry is spontaneously broken and are plentiful in string com- 
pactifications. The idea of using a PNGB as inflaton was first put forward in ||; in this 
minimal realization a single axion is present; it turns out that this model can produce a 
sufficiently flat spectrum of perturbations only if the axion decay constant / is above the 
Planck scale. This regime may be impossible to obtain in a controlled way: one can expect 
that gravity breaks the axionic shift symmetry at a smaller scale and it has also been con- 
jectured that / > M p cannot be realized in string theory. These problems have been solved 
by a number of controlled realizations in which an axion, or a combination of axions, with 
an axion scale / a few orders of magnitude smaller than M p , behaves effectively as a large 
field inflaton |, g^, 0, H U, 0, U • 

Axions are coupled to gauge fields through jipFF. More generally, from the perspec- 
tive of effective field theory, a coupling (pFF must be included whenever (p is pseudo-scalar. 
The dimensionless coupling a is a model-dependent quantity, however, from an field the- 
oretical point of view, one does not expect it to be <C 1. In this work we studied the 
phenomenological signatures induced by this coupling, keeping a/ f as a free parameter. 
The crucial point of our work is that, in presence of this coupling, the motion of (p induces 
a tachyonic growth for one polarization 5 A of the gauge field. This amplification is most 
important for modes with physical wavelength comparable to the horizon A ~ l/H. These 
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produced gauge quanta inverse decay to produce inflaton fluctuations of a comparable 
wavelength, consistent with the general expectation from causality/locality. These inflaton 
modes then leave the horizon, and the corresponding density fluctuations, Cinv.decay) become 
frozen and contribute to the observable cosmological perturbations. An analogous process 
leads to production of gravity waves. These scalar and tensor perturbations add up inco- 
herently with those generated by the expansion of the universe (the standard "vacuum" 
modes, ( vac and h va , c ). 

We find that the amplitude of the perturbations generated by the inverse decay is 
an exponentially growing function of a/f. These modes dominate over the vacuum ones 
for a/f > 10~ 2 Mp 1 (the precise value depending on the inflaton potential). Due to the 
exponential sensitivity on the coupling, there is only a small range in a/f for which both 
Cinv.decay; and Cvac are relevant. For smaller values of a/f this new effect is completely 
negligible, and the standard results are recovered. For larger values, drastically new pre- 
dictions are obtained. In particular, the main characteristic of Cinv.decay is that they are 
highly nongaussian: this is due to the fact that two gauge quanta participate in the inverse 
decay, and that the initial distribution of these quanta is itself gaussian (loosely speaking, 
Cinv.decay behaves as the square of a gaussian field, which is obviously not gaussian). As 
a consequence of this effect, the nongaussianity parameter f^L also grows exponentially 
with a/f in the region of parameter space for which both Cvac and Cinv.decay are comparable. 
When Cvac can be neglected, /tvl saturates to about 8, 600 in the equilateral configuration, 
well beyond the current WMAP limit f^j^ < 266. We find that the WMAP limit allows 
only values of a/f for which Cinv.decay can contribute < 10% to the power spectrum (see 
Figure ^|). In this regime, also the spectrum of gravity waves produced by the inverse decay 
is much smaller than that from vacuum. 

We have seen that nongaussianity from axion inflation is greatest for the equilateral 
configuration. This is related to causality /locality of the underlying inverse decay produc- 
tion mechanism. To understand this, intuitively, recall that a mode Cfc is sourced by the 
inverse decay of gauge perturbations of comparable wavelength. Consequently, correlation 

(CfciCfc 3 Cfe 3 ) K J d 3 qid 3 q 2 d 3 q 3 {5A qi 5A kl ^ qi 5A q2 5A k2 ^ q2 5A q3 5A k3 _ qa ) (7.1) 

is suppressed between modes of very different scale. We found that the shape of nongaus- 
sianity produced by this mechanism is very well reproduced by the equilateral template 
that has been widely studied in the literature. Quantitatively, the "cosine" between the 
two shapes (a measure of how much the two shapes coincide [^7]]) is = 0.93. Therefore, it 
is sensible to use the limits on f^ 11 (obtained with the use of the equilateral template) 
to probe the current mechanism. On the other hand, the specific shape that we have 
computed is nevertheless distinct from the equilateral template and this fact might be 
useful to distinguish axion inflation from other models (for instance, the shape we found 
has a significant overlapping also with the orthogonal template). This issue can be more 
thoroughly explored when / if a nonzero value for Z^" 1 will be found in the data. The 
nongaussian signature of axion inflation, together with the requirement of standard results 
for the spectral tilt and the tensor-to-scalar ratio tensor of large field inflation, will allow 
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to falsify the mechanism (or, at least, the simplest version that we have studied here) in 
the near future. This is due to the fact that large field inflationary model provide much 
larger - and detectable - values of n s — 1 and of r than many other classes of models. 

We stress that the nongaussianity which we have studied is very different from the 
so called resonant nongaussianity [^7|, 28], that has been discussed previously [|^, 59 



|60| , 44] in the context of axion monodromy inflation [12]. Axion monodromy is one of the 
particle physics realizations of axion inflation; in this model the periodic potential typical of 
axions is only a subdominant term in the inflaton potential, and it provides a nongaussian 
modulation of the inflaton perturbations. Depending on model parameters, either this 
effect or inverse decay processes may dominate the bispectrum. 

It is remarkable that, in the mechanism we have studied, the large nongaussianity 
is obtained in a rather minimal way (simply by considering a coupling of the pseudo- 
scalar allowed by the symmetries of the model, and therefore expected in an effective field 
theory context). Nongaussianity is a measure of the strength of interactions of the inflaton, 
which are typically constrained to be small by the requirement of a flat potential. To 
obtain observable nongaussianity, previous studies have invoked non-standard field theories 
(involving small sounds speed [23] or higher derivatives ||), or initial conditions [23, 24, 
|25| . |26[ , potentials with sharp features [p^, [28]], dissipative effects ||, fine-tuned inflationary 
trajectories (29] or post-inflationary effects (such as preheating [30, |3l|). The mechanism 
we have studies circumvents the common lore result in a very novel way: the interaction 
that gives rise to large nongaussianity, eqn. Ql.3j ), does not play any role in the background 
dynamics and is thus unconstrained by the requirement of slow roll. At the same time, 
the effect of this interaction (namely, the amplification of gauge quanta) persists during 
the entire inflationary phase, leading to a (nearly) scale invariant signature (as opposed to 
sporadic episodes of particle productions |3^, |33|, [34]] that would lead to a highly localized 
nongaussianity in momentum space). 

Finally, it is also remarkable that current observational limit on nongaussianity already 
place a surprisingly stringent bound on pseudo-scalar interactions of the form ( |1.3|) . As 
a comparison, the bound on the coupling of the much lighter QCD (or QCD-like) axion 
to photons is a p hotons// ~ O (l0~ n ) GeV -1 , from energy loss in stars J74|. On the other 
hand, the mechanism that we have studied provides the bound a/f^O (lO -16 ) GeV -1 
on the coupling of the pseudo-scalar inflaton with any gauge field. This provides a unique 
window for constraining - or perhaps probing - a large class of inflationary models. As we 
have reviewed in section ^, there are a variety of interesting multi-field models involving 
axions, many of which can be realized in string theory. While most of these studies provide 
the value of axion decay constant / and of the potential V((p), comparatively less attention 
has been devoted to compute the dimensionless coupling, a. We believe that our findings 
provide a strong motivation to undertake such a study. 
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APPENDIX A: Gauge Field Mode Functions 



In this appendix we discuss the solutions of equation ( 2.1C| ), describing the unstable growth 
of gauge field fluctuations in the background of the slowly rolling inflaton (hence, £ can be 
taken constant at leading order in the slow roll parameters). We require that the gauge 
field is initially in the adiabatic vacuum: A±(r,k) = e~ lkr j\[2~k for kr — > — oo. The 
mathematical properties used here can be found in chapter 14 of [fT5| . The solution of 
(2.10) which satisfies this condition may be expressed in terms of Coulomb functions: 

A+(t, k) = -i= [G (£, -kr) + iF & -kr)] (A-l) 



2k 

The production of gauge fluctuations is only interesting in the region of phase space —kr-^i 
2£ and when e*^ 3> 1 (sere the discussion in subsection 2.2). In this regime equation ( |A-1| ) 
may be very well approximated in terms of the modified Bessel function of the second kind, 
K u (z), as 



A+(r,k) 



-It 



7T 



e^Ki \2y/-2£kr 



(A-2) 



The growth of the modes (A-2) saturates deep in the IR: for —kr we have A + — > 
e 7r ^/(2 v / 7rA;£) so that the physical electric and magnetic field vectors ( |2.7| ) decay suffi- 
ciently far outside the horizon. An inspection of the solutions shows that the interesting 
physical effects (for instance, the production of £ 2 and £ 3 correlators) take place in the 
region (8£) _1 <, — kr < 2£ of phase space. In this regime we can take the large argument 
asymptotics of the Bessel function in ( [A- 2 ) to obtain a very simple representation of the 
modes: 



A + (r,k) 



A'(r,k) 



'2k 



2£ J 



2fc£ 



A+(r,k) 



(A-3) 



Throughout the majority of this paper we employ the representation ( |A-3| ) of the 
modes, for brevity of exposition. However, we have verified that none of our results changes 
significantly if we use the more accurate expression ( |A-2 ). Formally, the only effect of using 
(| A- 2D rather than (|A-3| ) is that, for any two (rescaled) momenta qi and q2, the quantity 



« I 1 / 2 _1_ In I 1 / 2 

qi| + I<i2| 



dx (sinx - xcosx) e - c (l qi l 1/2+ l q2 l 1/2 )v /i (A-4) 

Jo 



entering in eqs. ( |3T7| ) and ( |3^9| ) gets replaced by 



z ^ _ 2c |q 1 | 1 / 4 |q 2 | 1 / 4 
7T (qi) 1 / 2 + | q2 1 1//2 



dx 



x \smx — xcosx 



qi 



c |qix| 



l ' 2 ) K, 



c|q 2 x 



1/2 



(A-5) 

+ |q 2 | 1/2 ^i (cl^xl 1 / 2 ) K fclqixl 1 / 2 



We verified that the results obtained with this replacement are in excellent agreement 
with those presented in the main text. Specifically, for the nongaussianity shape we ob- 
tained identical values to those presented in Table 1 for the nongaussianity shape; for the 
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nongaussianity parameter /j^l" 1 we obtained values consistent within a few percent with 
those shown in Figure ||. This explicitly confirms that the approximated expressions (|A-3j ) 
are adequate for our analysis. 

APPENDIX B: Perturbing the Klein-Gordon Equation 



In section |5| we have derived a closed- form expression for Klein-Gordon equation ( 5.15| ) 
for the inflaton fluctuation 5ip in the theory (fOI). This was derived by employing the 



ADM decomposition of the metric (5.3) and integrating out the lapse N and shift N % 
functions. However, one could arrive at precisely the same result by working directly with 
the equations of motion. Following [EE], we expand the metric up to second order in 
perturbation theory as 



a 2 (l + 20i+20 2 ) 



a 2 d i (B 1 +B 2 ) 



goo 

9oi 

9ij — " u tj 

Similar, we expand the scalar and gauge fields as 



a 2 5j 



(B-l) 
(B-2) 
(B-3) 



<p(t,x) = 4>(t) + 6i<p(t,x) + 5 2 tp(t,x) 
A^(t,x) = (0,5 1 A i (t,x)+5 2 A i (t,x)) 



(B-4) 
(B-5) 



In this gauge, the curvature perturbation on uniform density hypersurfaces is £ 



— £[5<p. The equation of motion for the scalar field is 



(B-6) 



We expand (B-6) up to second order in perturbation theory, using the Einstein constraint 
equations to eliminate the metric fluctuations in order to close the system. At linear order 
we find 



dl + 2Ud T - V 2 + a 2 m 2 



M 2 
p 



5 x ip = 



to leading order in slow roll. At second order in perturbation theory we find 



(B-7) 



d 2 d 

di + m d: 



V 2 + a 2 m 2 



a'2 



d 2 tp 



a 1 



fa 



2 e ijk 5iA' i d j (5 1 A k ) 



+ 



2a 2 UMf } 



+ 0[(5 W ) 2 ] (B-J 



to leading order in slow roll. These results coincide with those reported in section [5], 
which were derived starting from the action of the perturbations. Besides being a check 
on our algebra, this agreement strengthens the conclusions of helping to establish the 
consistency of the two most popular approaches to nonlinear cosmological perturbation 
theory. 
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